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ICCUm TV  CLASSIFICATION  OF  THIS  FAOe  fWfc*a  0«f«  Emtm++G 


Bounds  on  Distributions  Arising  in  Order  Restricted  Inference:  The 
Partially  Ordered  Case 

ABSTRACT 

In  testing  hypotheses  involving  order  restrictions  on  a  collection 
of  parameters,  distributions  arise  which  are  mixtures  of  chi-squared  or 
beta  distributions.  In  general,  the  mixing  coefficients  are  quite 
intractable  even  for  a  moderate  number  of  populations.  Stochastic  upper 
and  lower  bounds  are  obtained  for  mixtures  which  arise  in  Bartholomew's  tests 
for  homogeneity  of  normal  means  with  the  alternative  restricted  by  a  quasi 
ordering.  These  bounds  are  applicable  in  the  dual-testing  situation,  that 
is  in  testing  the  order  restriction  as  a  null  hypothesis;  in  testing  order 
restrictions  in  exponential  families,  Poisson  intensities  and  multinomial 

ru^ 

parameters;  and  in  some  nonparametric  settings.  They  can  also  be  applied 
to  obtain  the  least  favorable  configuration  for  testing  equality  of  two 
multinomial  populations  with  a  stochastic  ordering  alternative. 


‘  '55 i on  For 
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1  INTRODUCTION.  Suppose  one  is  interested  in  testing  hypotheses  about  a 


vector  of  normal  means,  p  =  (  p  1  ,  p^ »  •  •  •  •  ) »  and  it  is  known  apriori  that 

these  means  satisfy  certain  restrictions,  such  as  ^  i  or  ij.  1 

0  for  i  =  1,2 . k.  Likelihood  ratio  tests  (LRTs)  which  make  use  of  such 

order  restrictions  have  been  derived  under  various  assumptions  on  the 
covariance  structure. 

Bartholomew  (1959,  1961)  considered  testing  homogeneity  subject  to  the 
restriction  that  the  p^  satisfy  a  partial  ordering.  Let  <<  be  a  quasi  order 

on  T  =  {1,2,...,kJ,  that  is  <<  is  reflexive,  ie.  i  <<  i  for  all  i  e  f,  and 
transitive,  ie.  i  <<  j  <<  £  implies  i  <<  £.  (A  partial  order  is  also 
antisymmetric,  ie.  i  <<  j  <<  i  implies  i  =  j.)  If  u.  <  u.  for  all  i  <<  j , 


then  p  is  isotonic  with  respect  to  (wrt)  <<•  Since  k  =  1  is  not 
interesting,  we  suppose  k  12.  For  j  =  1,2,...,n  and  all  i  e  r,  let 

2  2  2 

( A 1 )  X .  .  be  i ndependent ,  X..  ~  N(p.,o.)  and  o .  =  a . o_  with  a .  known  and 
ij  ij  ii  l  l  0  l 

Op  =  1  if  known. 


Bartholomew  studied  the  LRT  of  Hq:  p1  =  =  —  =  uk  versus  Hj-Hq  with  H1  : 

p  is  isotonic.  He  showed  that  the  null  distribution  is  a  mixture  of  chi- 

2 

squared  or  beta  distributions  depending  on  whether  oQ  is  known  or  not. 

Under  a  more  general  assumption  about  the  covariance  structure,  Perlman 
(  1  969)  considered  the  LRT  of  versus  H^:  p  is  not  isotonic  and  obtained 

bounds  for  the  significance  level  of  the  test.  Assuming  (A1),  Robertson  and 
Wegman  (1976)  proved  that  is  least  favorable  for  testing  H ^  versus  and 

showed  that  the  distribution  of  the  LPT  statistic,  under  H ,  is  again  a 
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mixture  of  chi-squared  or  beta  distributions  with  the  same  mixing 
coefficients  as  in  Bartholomew's  test. 

The  mixing  coefficients  depend  upon  the  quasi  order,  <<,  and  w  = 
(w  ,w9, .w^)  where  w.^  =  n./a^  ^or  *  c  r*  Even  in  the  case  in  which  <<  is 

a  total  order  (ie.  i  <<  j  or  j  <<  i  for  all  i,j  e  D,  these  coefficients  are 
quite  intractable  if  the  weights,  w.,  are  not  equal  and  k  >  5,  cf.  Barlow  et 

al.  (1972,  p-  1*12).  (The  equal-weights  mixing  coefficients  for  a  total 
order  are  given  in  their  Table  A. 5.)  For  the  totally  ordered  case  with 


unequal  weights  and  oQ  known  or  unknown,  the  distributions,  under  H^,  of  the 

LRT  statistics  for  HQ  versus  H1  and  H1  versus  have  been  studied. 

Approximations  were  developed  by  Chase  (197*0,  Syskind  (1976)  and  Robertson 
and  Wright  (1983).  Upper  and  lower  bounds  for  their  significance  levels  are 
given  in  Robertson  and  Wright  (1982). 

For  the  partially  ordered  case,  less  is  known  about  the  null 
distributions  for  these  tests.  An  iterative  scheme  for  computing  the  mixing 
coefficients  is  discussed  in  Barlow  et  al.  (1972,  p.  139).  They  also  give 
exact  expressions  for  these  coefficients  for  arbitrary  weights,  k  <  and 
the  following  partial  orders: 

(1)  the  simple  tree;  1  <<  i  for  i  =  2,3,-..,k  (the  equal-weights 
coefficients  for  a  simple  tree  are  given  in  their  Table  A. 6)  and 

(2)  the  simple  loop;  1  <<  i  <<  k  for  i  =  2,3,...,k-1. 

There  are  other  partial  orders  that  are  of  practical  importance,  too.  For 
instance,  if  is  the  mean  response  of  a  dependent  variable  with  one 

independent  variable  at  level  i  and  the  other  at  level  j  and  u. .  is 

i  J 

nondecrea s i ng  in  i  with  j  held  fixed  and  vice  versa,  then  u  =  (u  )  is 


isotonic  wrt  <<,  the  coord inatewise  partial  ordering.  (This  example  is 

discussed  in  Brunk,  Ewing  and  Utz  (1957)  and  Hanson,  Pledger  and  Wright 

(1973).)  The  unimodal  ordering,  p,  i  p„  S S  p,_  2  p.  ,  Z  —  £  p.  with  1 

1  2  h  h+ 1  k 

<  h  <  k  is  also  of  interest.  There  is  much  to  be  learned  about  the  mixing 
coefficients  for  such  partial  orders.  Robertson  and  Wright  (1985)  and 
Wright  and  Tran  (1985)  give  approximations  and  bounds  for  the  case  of  a 
simple  tree. 

In  Section  2,  we  show  how  the  quasi  ordered  case  is  reduced  to  the 
partially  ordered  case  and  obtain  sharp  bounds  for  the  significance  levels 
of  the  LRTs  of  HQ  versus  and  H1  versus  which  are  independent  of  the 

weights,  w^ .  In  Section  3,  some  applications  of  these  bounds  are  discussed. 

Kudo  (1963)  considered  a  closely  related  problem,  that  of  testing  H^:  6 

=  0  versus  HJ-H^  with  H| :  6  £  0,  where  9  1  0  means  each  coordinate  of  6  is 

nonnegative.  He  derived  the  LRT  based  on  a  .random  sample  from  a 
multivariate  normal  distribution,  ie. 

(A2)  Y.  -  ( Y Y _,..., Y .  )  are  i.i.d.  N(6,E)  for  i  =  1,2,...,n  with 

l  1 1  1 2  i  p 

I  nonsingular, 

2 

for  the  case  of  known  E.--T«eusch  (1966)  treated  this  problem  with  l  =  oQ/\ 

and  A  known.  (See  the  comments  in  Section  12  of  Shorack  (  1967).)  Perlman 
(1969)  considered  the  case  of  unknown,  nonsingular  £.  In  fact,  Perlman 
derived  the  LRTs  for  0=0  versus  0  r  S-{0)  and  0  c  C  versus  9  f  C  with  C  a 

positively  homogeneous,  one-sided  subset  of  (see  that  paper  for 
definitions)  and  obtained  sharp  bounds  for  the  significance  levels  of  these 


tests.  Some  of  the  testing  problems  discussed  earlier  can  be  placed  in  this 


from  -work . 


for-  instance,  for-  the  simple  tree  with  variances  known,  set  n 


1,  p  =  k-1  ,  Y.  .  -  X.  +  1  ~  xi  for  J  “  1,2,.  ...p  and  C  the  p-d  imens  iona  1 
J  J 

positive  orthant.  The  total  order  and  unimodal  order  can  be  treated 
similarly.  However,  if  the  variances  are  unknown  and  the  sample  sizes  are 
not  the  same,  or  if  this  differencing  yields  a  nonsingular  E  (see  the 
discussion  of  the  simple  loop  in  Barlow  et  al.  (1972,  p.  176)),  then  this 
technique  does  not  seem  to  work.  On  the  other  hand,  (A 2)  includes 
situations  not  included  in  (Al).  The  bounds  given  here  are,  for  some 
partial  orders,  tighter  than  the  bounds  obtained  from  Perlman's  work,  but 
that  is  because  they  depend  on  the  special  covariance  structure  implied  by 
(Al).  Perlman  shows  that  for  any  C,  which  contains  a  p-d  imensional  open 
set,  his  bounds  are  sharp. 

2  TAIL  PROBABILITY  BOUNDS.  Throughout  this  section  we  assume  that  (Al) 

2 

holds  unless  stated  otherwise.  If  cQ  is  known  the  LRT  of  HQ  versus  H1 
rejects  for  large  values  of 


with  u 


). k  .  w  .  X  .  / 1 k  w. 

i=1  l  l  i=1  l 


and  u 


E  ( X | H 1 ) ,  the  projection  of  X  = 


(X  ,X^,...,X  )  onto  H1  with  respect  to  (wrt)  the  distance  determined  by  the 

rk 

inner  product  (x,y)  =  ).  , w.x.y..  Furthermore,  under  H. 

w  l  =  1  l  l  l  0 

<?•»>  PwUq,  5  C]  =  );^1PU.k;w)P[X'-_1  >  cl 

')  _ 

with  Xn  5  0  and  P(5.,k;w)  the  probability,  under  H  ,  that  E  ( X  |  H  )  has 
u  'j  w  i 


exactly  l  distinct  values  (cf.  Barlow  et  al.  (197 2,  Chapter  3))-  Of  course, 


6 


XQ1  can  be  written  as  ||P-w||  .  in  which  p  also  denotes  a  k-d  imensional 

vector  with  each  component  equal  to  p.  Barlow  et  al.  also  discusses  the  LRT 

2 

of  Hq  versus  if  oQ  is  unknown.  It  rejects  HQ  for  large  values  of 


3.  • 


Under  H„, 
0 


(2.2)  PIZ  >  o]  ■  2  o) 

-U-1 

with  B  a  beta  variable  with  parameters  a  and  b  (B.  .  =  0)  and  N  » 
a ,  b  0 ,  b 

rk 

l .  ,  n .  . 

l  =  1  l 

2 

If  o  is  known,  the  LRT  of  H  versus  H  :  p  is  not  isotonic  rejects  for 


large  values  of 


hi  ■  n-,  VW"  ■ 


is  least  favorable  within  H1 ,  and  under  HQl 


(?.3)  Pwl'xf2  ^  c]  =  Ij  =  1P(i,k;w)P[x^_t  >  cj. 


If  o.  is  unknown,  the  LRT  of  H  versus  H  rejects  for  large  values  of 


r  ,  w.(X.-p.)'V|)k  ,  w  .  ( X  .  -  p  .  Y  *  Q] 

1=1  1  L  1  1=1  1  1  1 


,  ,  n  . 

,  k  -  1  r  i 


with  Q  =  a.  J . =1  (X.j-X.  )  .  HQ  is  again  least  favorable  within  1^  and 


under  Hq, 


7 


(2.D  P  LE2  2  c)  -  £  ,PU.k;u)PtB,  ,  2  o]. 

“  12  *■'  |(k-«).i(»-k) 

(See  Robertson  and  Wegman  (1978).) 

We  seek  upper  and  lower  bounds  for  the  significance  levels  given  by 

(2. 1) -(2.4).  The  techniques  from  Perlman  (  1969)  and  Robertson  and  Wright 
(1982)  will  be  combined.  In  the  latter  reference,  it  was  shown  that  for  the 
total  order,  1  <<  2  <<  ...  <<  k,  an  upper  (lower)  bound  for  (2.1)  and  (2.2) 

k-i  -k+1 

((2-3)  and  ( 2 .  t| ) )  is  obtained  by  replacing  P(J.,k;w)  by  A(2.,k)  =  (  )2  , 

1  <  i  <  k,  in  those  expressions,  and  a  lower  (upper)  bound  for  (2.1)  and 

(2.2)  ((2.3)  and  (2.U))  is  obtained  by  replacing  P(£,,k;w)  by  B(£,,k)  = 

|  I(1  9}(g.)  with  I  ( 2. )  =  1  if  l  e  A  and  IflU)  =  0  if  l  j  A.  They  also 

established  the  sharpness  of  these  bounds.  The  bounds  resulting  from  the 
13(8., k)  were  obtained  by  Perlman  (1969)  and  the  mixing  coefficients,  A(£,k), 
were  encountered  by  Dykstra  and  Robertson  (1982,  1983)  in  related  testing 
problems . 

If  Y  is  a  binomial  variable  with  parameters  k-1  and  ^ ,  then  the 

probability  function  of  Y  +  1  is  A(£,k),  and  so  we  refer  to  the  {A(£.,k)J  as 
shifted  binomial  probabilities.  Because  the  bounds  of  interest  here  involve 
convolutions  of  shifted  binomial  probabilities,  we  give  the  following 

k-b  — k ♦ b 

ciefinition.  For  positive  integers  b  <  k,  set  A^(2,k)  =  (  ^)2  for  £  = 

1,2 . k.  (Of  course,  A^(Z,k)  =  0  for  ?.  <  b  and  A(fc,k)~  Aj(£,k).) 

The  LPT  statistics  and  their  null  distributions  given  in  this  section 

were  derived  for  partial  orders.  However,  if  u  is  understood  to  be  the 

projection  of  X  onto  the  cone  corresponding  to  a  quasi  order,  <<,  then  the 
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arguments  given  in  the  references  above  show  that  the  LR  procedure  gives 
rise  to  the  same  test  statistics,  Xq1»X^2,E01  and  ^12'  Dykstra  (1981) 

discusses  the  computation  of  u  for  an  arbitrary  quasi  order.  The  argument 
given  to  show  that  HQ  is  least  favorable  within  in  the  partially  ordered 

case  is  also  applicable  for  an  arbitrary  quasi  order.  In  the  next  remark, 
we  show  that  (2.1)-(2.iO  are  also  valid  for  quasi  orders.  The  proof 
consists  of  reducing  <<  to  a  partial  order  and  applying  the  known  results. 
This  reduction  also  shows  that  we  only  need  to  obtain  bounds  for  the 
partially  ordered  case. 

If  <<  is  a  quasi  order,  then  i  =  j  if  and  only  if  i  <<  j  <<  i  is  an 

equivalence  relation  and  induces  a  natural  partial  order,  <<*,  on  the 

equivalence  classes  {e  :  a  =  1,2,  —  ,k*}.  Vie  identify  a  with  e  ,  think  of 

a  a 

k  * 

<<*  as  a  partial  order  on  {1,2,... ,k*}  and  let  H*  denote  the  vectors  in  R 

isotonic  wrt  <<*.  For  1  <  a  S  k#,  set  w*  =  7.  w.  and  X*  =  ).  X./w*; 

a  lee  l  a  ice  i  a 

a  a 

let  ij*  denote  E  „(X*|H*);  and  let  P(S.,k*;w*)  denote  the  probability  under 

homogeneity,  that  u*  contains  exactly  1  distinct  values.  (Of  course 
P(". ,k*;w*)  =  0  for  k*  <  9,  S  k.) 

IP-mark.  Under  homogeneity,  (2.1 )-(2.U)  are  valid  for  quasi  orders  and  their 
right  hand  sides  are  unchanged  if  P(<.,k;w)  is  replaced  by  P(  8. ,  k*  ;w* ) . 

Proof.  Clearly  u*  =  >  ,  w*X*/)  '  w*  =  q.  Since  E  (X  H,  )  is  constant  on 
-  ^  Q  n  ‘-a«  l  a  w  1  1 

each  e  ,  it  must  agree  with  E  _  ( X  *  I H  * )  o,n  these  equivalence  classes.  Hence, 
a  w  *  1 


P(£,k;w)  *  P(i.,k*;w*)  for  t  ^  k,  and  consequently,  the  second  conclusion  is 
established.  Next,  we  observe  that 


(?.5)  Xq -|  3["  =  1  I*!,  lice  w.(w*  - 

a 


rk* 


*i2  Yk* 

J  >  =  Cl 


w*(u*  -  u*)  -  U 
a  a 


and  since  u  is  constant  on  e  , 

a 


(2.6)  xfj,  Hlj  =  1  w.(X.  -  uJ2  =  Q'  +  V 


?  ~k#  —  —  2 

with  Q'  —  7  ,  I.  w.(X.  -  X*)  and  V  ■=  £  ,  w*(X*  -  u*)  .  We  assume  H 

—  a=  1  ‘'ice  li  a  a=1  a  a  a 

a 


0 


for  the  remainder  of  the  proof.  It  follows  from  (2.5)  and  Bartholomew's 


work  that  (2.1)  holds  for  quasi  orders.  Clearly,  C,  C'  and  X*  are 

2 

independent  and  Q'  “  x  (k-k*).  Conditioning  on  Q'  and  applying  Corollary 
(ii.2)  of  Robertson  and  Wegman  (1978),  we  see  that 

Pw[^12  "  c]  =  Cl  PCt.k^.-wMPtx2.,,  *  c] 


— ?  —2  -2 

and  so  (2.3)  also  holds  for  quasi  orders.  Rewriting  EQ1  as  X01/(Q  *  XQ1  * 

_2  -2  —2  2 

X 1  p )  ,  recalling  that  Xq1  =  U  and  x12  =  V+Q'  and  noting  that  Q  +  Q'  x  (N- 

k*),  wc  see  that  E^  has  the  same  distribution  as  the  correspond i ng  test 
statistic  for*  th:'  partial  order,  <<* ,  with  a  total  simple  sir..*  of  N.  Thus, 

_2 

applying  Bartholomew's  results  on  E  ,  we  see  that  (2.71  holds  for  an 

_ 'I 

arbitrary  quasi  order.  dinm*  0  1  i  1,  wo  consider  e  r  (0,1)  and  note 


that 


10 


P[E2]2  >  c]  =  PCxfo/Q  *  C/(1-C)]  =  P[v  >  cQ/(  1-c)  -  Q ’  ] . 

Applying  the  result  of  Robertson  and  Wegman  for  V  to  the  above  probability, 
conditioned  on  Q  and  Q' ,  we  obtain 

2 

—2  rk*  xk*-£  + 

P[E  >  c]  =  [  PU.k«;w*)P[-5 -  >  c] 

^  v  +Q'+0 

Xk*- Z  V  V 

2 

with  xk*  Q  and  Q'  independent.  Thus,  (2.4)  holds  for  quasi  orders,  too. 
The  proof  is  completed. 


-2  -2 

Because  the  null  distribution  of  Xqj  and  fo r  a  quasi  order,  <<,  is 

the  same  as  for  the  associated  partial  order,  <<*,  bounds  for  a  partial 
order  are  also  bounds  for  any  quasi  order  that  reduces  to  the  given  partial 
order.  Considering  the  proof  of  the  Remark,  it  is  clear  that  if 

v  k  *  2 

aJP[x(<)(_i  2  c]  provides  a  bound,  either  upper  or  lower,  for  the  tail 
-2 

probability  of  x12  under  HQ  for  the  partial  order  <<*,  then 


vk* 


dtPl-xk-t 


o]  (1 


ck* 

i=1 


VLB1  1 

^(k-£) ,^(N-k) 


c]) 


is  in  analogous  bound  for  the  tail  probabilities  of  x1  >(1-^)  under  H  for 

the  quasi  order  <<.  Hence,  we  need  only  obtain  bounds  for  the  partially 

eref  rro  c  :  :;e  . 

W< •  suppose  that  <<  is  a  partial  order  and  consider  upper  (lower)  bounds 
for  (2.1)  and  (2.2)  ((2.3)  .and  ( 2 .  <1 ) ) .  First,  we  need  the  following. 

De  f  l  n  i  t.  i  on:; .  Two  distinct  elements  i  ,  j  c  f  are  comparable  if  i  <<  j  or  j  << 
i.  Otherwise,  they  ere  noncomparable. 


A  subset  of  I'  is  a  chain  (antichain)  if  each  pair  of  distinct  elements 


in  the  subset  is  comparable  (noncomparable). 

The  breadth,  b,  of  r  is  defined  to  be  the  maximal  cardinality  of  any 
antichain  in  I',  ie.  b  =  max{card .  (A) :  A  an  antichain  in  r}. 

The  following  result  is  given  in  Crawley  and  Dilworth  (1973,  p.  3)- 
Theorem  2.1.  Let  (r,<<)  be  a  partially  ordered  set  with  breadth  b.  T  can 
be  written  as  the  disjoint  union  of  b  sets  which  are  chains  wrt  <<. 

Theorem  2.2.  Let  (r,<<)  be  a  partially  ordered  set  with  breadth  b,  and  let 
u  t:  11^ .  For  any  w  with  w.  >  0  for  each  i  c  f, 


PWC Tq,  i  c]  S  TI,,«b(t.k)P[x^1  2  c]. 


(2.7b)  P  [E.3  >  c]  <  lk.  A  (£.,k)P[B*  .  >  c], 

w  01  11=1  b  lu-i).l(N-U 


(2.70  Pulxfp  a  O]  >  lj=1  Ab(Z,k)P[x^_s  2  c]  and 


(2.7d)  P  [k'  2  c]  >  lk  A  ( l  ,k  )P[B  .  2  c] 

i(k-!l),i(li-k) 


Proof.  Usinft  (7.6)  and  (7.7)  of  Barlow  et  nl.  and  the  fact  that  Hq*  it 


’  '*  ’';i i  '  7  :>hown  that.  1  I  X  -  n  I  I  w  =  I  !  Kw  ( X  |  i5 ,  )  -  ullf  +  I  !  X  -  Ew(x|H,)|lw- 
Furthermore ,  if  HJ  correspond:’,  to  some  quasi  order  on  f  and  H  j  ,  then 


| '  - 


|  X  -  K  ( X | H ! ) |  |  and  consequent 1 y  f  j  | K  t ( X | H  ) 


I  I  F  (X  I  It;  )  -  „!  T  .  L«'t.  F  ,  F  ,,...,  r  be  disjoint  chains  with  b  ,  T  =  f;  l«t 

1  1  w  1  1  1  w  1  1)  a=  1  a 

>:  ■  ’  be  defined  by  i  <<’  j  provided  i,j  >:  ('  for  some  a  and  i  <<  j;  and  let 


k;  the  other  is  similar.  Since  k  is  an  interior  element  there  exists  c  ^  k 
with  k  <<  c.  Of  course,  c  e  A  or  c  e  B,  but  either  case  contradicts  A  +  B  = 

r  -  (k). 

It  is  easy  to  show  that  r  -  { k )  has  i-1  interior  and  e  exterior 

(k-1  ) 

elements.  By  the  inductive  hypothesis  P(S.,k-1;w  )  can  be  made  close  to 

A(5.,e).  The  desired  conclusion  follows  readily,  and  the  proof  is  completed. 

3 _ A P Pkl^ATHlNS .  The  bounds  given  in  Section  2  provide  conservative 

critical  values  or  p-values  for  testing  versus  H1  and  versus  with 

p  1  , p ^  .  •  •  •  .  u ^  normal  means.  (Some  particular  partial  orders  will  be 

discussed  in  this  section.)  Robertson  and  Wegman  (1978)  show  that  the  chi- 
bar-squared  distributions  arise  as  large  sample  approximations  for  the  LRTs 
of  trends  in  parameters  in  an  exponential  family.  Robertson  (1978) 
considers  testing  for  trends  in  multinomial  parameters.  Also,  nonparametr ic 
tests  for  trends  can  bo  obtained  by  replacing  each  observation  by  its 
ranking  among  all  the  observations  (cf.  Shirtey  (  1977)).  Chi-bar-squared 
distributions  arise  in  this  setting,  too.  The  bounds  established  here  can 
be  applied  in  each  of  these  three  cases. 

For  a  total  order  and  the  simple  tree,  such  bounds  were  useful  in 

'h>t*:rrr  i  n  ing  approximations  for  the  x*  ■'snd  K*  d i otr ibut ions ,  in  studying  the 

_ j  _ /> 

robustness  of  the  equal-weights  (w  =  w .  =  ...  «  w  )  x'  and  K‘  distributions 

and  for  identifying  weight  sots  for  which  the  equa  1 -we ights  distributions 
might  not  provide  very  accurate  approximations.  It  is  anticipated  that  the 
bounds  given  here  will  be  used  similarly.  For  instance,  in  the  case  of  a 
simple  loop,  we  see  from  the  proof  of  Theorem  2.*J  that  if  w1  and  w^  are 


2*4 


exterior  elements,  which  is  at  least  two  since  r  is  finite  with  at  least  two 
elements.  If  T  is  indecomposable  with  no  interior  elements  and  e  =  2,  then 

P ( 1  , 2 ; w  )  =  P (2 ,2 ; w )  =  1. 

Hence,  we  suppose  that  the  conclusion  holds  for  indecomposable  sets 
with  no  interior  elements  and  e-1  exterior  elements.  Consider  r,  an 
indecomposable  partially  ordered  set  with  no  interior  elements  and  e 
exterior  elements.  By  Lemma  2.7  there  exists  j  e  f  with  f  -  { j  } 

indecomposable,  and  it  is  easy  to  see  that  r  -  { j }  has  no  interior  elements 
and  e-1  exterior  elements.  For  convenience  relabel  so  that  j  =  k.  Because 
T  is  indecomposable,  j  is  not  isolated,  and  by  Lemma  2.3,  we  can,  by  making 
w^  small,  make  {P(2.,k;w)}  uniformly  close  to  the  convolution  of  {P(JL,k  — 

1;w^k  '))}  with  I  j Q  1  j ( £ ) )  -  By  the  inductive  hypothesis  (P(£,k- 

1,-w^k  '))}  can  be  made  uniformly  close  to  {A(£.,c--1)l  and  the  desired 
conclusion  follows. 

Returning  to  the  primary  induction,  we  suppose  that  the  desired 
conclusion  holds  for  indecomposable  sets  with  i-1  interior  elements  and  that 
f  is  indecomposable  with  i  interior  elements  and  e  exterior  elements  (of 
course,  c»i  =  k).  For  convenience  suppose  that  k  is  an  interior  element. 
By  Lemma  2.6,  we  can,  by  choosing  w^  sufficiently  small,  make  P(?,,k;w)  close 

(  k  -  1  ) 

to  P  (  H  ,  k  - 1  ;  w  ).  Next,  we  show  that  r  -  { k )  is  indecomposable.  If  not. 

there  exist  A  and  B  with  A  *  B  =  r  -  lk[.  There  must  be  elements  in  both  A 
and  B  which  are  comparable  with  k,  for  if  not,  say  B  contains  no  elements 

comparable  with  k,  then  (A(J  fk])  ♦  B  =  f.  This  is  a  contradiction. 
Suppose  a  e  A  and  b  c  B  with  a  and  b  both  comparable  with  k.  Now,  either  a 

<<  k  and  b  <<  k  or  k  <<  a  and  k  <<  b.  We  consider  the  ease  a  <<  k  and  b  << 


Proof .  The  proof  is  by  induction.  If  k  =  2,  the  conclusion  is  obvious. 

Consider  k>2.  If  r  -  { k }  is  decomposable,  then  there  exist  A  and  B  with  A  + 
B  =  T  -  { k }  .  We  will  show  that  A  U  { k }  is  indecomposable  (throughout  the 
proof  we  mean  indecomposable  wrt  the  appropriate  restriction  of  <<).  If  not 

there  exist  A*  and  B*  with  A'  +  B'  =  A  U{k}.  Without  loss  of  generality, 
we  assume  k  ^  A*.  Since  A'  C  A,  we  know  that  every  element  in  A*  is 
noncomparable  with  every  element  in  B.  Because  every  element  of  A'  is 

noncomparable  with  every  element  in  B' ,  A*  ♦  (B  U  B' )  =  T,  a  contradiction. 
By  the  inductive  hypothesis,  there  exists  j  e  A  (J  { k }  with  A  U  (k)  -  { j } 

indecomposable.  Next,  wc  show  that  V  -  {j}  is  indecomposable.  If  not  A"  + 
B"  =  T  -  (j}.  Considering  separately  the  cases  j  =  k  and  j  f  k,  we  see  that 
( T  -  { j  } )  O  (A  U  (k)  -  {j})  ^  <ji.  For  convenience,  let  A"  fl  (A  O  { k }  -  { j  }  ) 
f  <p.  Because  A  U  (k)  -  ( j  }  is  indecomposable,  A"  All  { k }  -  { j ) .  It  is 
not  difficult  to  show  that  B"  C  B.  Every  element  in  B,  and  consequently 
every  element  in  B" ,  is  noncomparable  with  j.  Also,  every  element  in  A"  is 

noncomparable  with  every  element  in  B" .  Hence,  (A "  \J  { j  }  )  +  B”  =  T,  which 
is  a  contradiction. 

Theorem  2.8.  If  (T,<<)  is  indecomposable  with  e  exterior  elements,  then 
there  exists  a  sequence  of  positive  weights  w(n)  =  (w, (n) ,w0(n) , . . . ,w  C  n ) ) 

l  t  k 

for  which 

lin  P(t,k;w(n))  =  A(> t,e)  for  l  =  1,2,...,k. 
n  *” 

ol’  course,  A(i,e)  =  0  for  d  >  e. 

Proof .  The  proof  is  an  induction  on  the  number  of  interior  elements.  We 
first  show  that  the  conclusion  is  valid  for  indecomposable  T  with  no 
interior  elements.  This  part  of  the  proof  is  an  induction  on  the  number  of 


This  produces  a  partition  in 


and  the  probability  that  this 


C. 


f  (2) 

Ik  ’ 


partition  occurs  converges,  as  ♦  0,  to  the  probability  that 

{C1  .C2*  '  *  * ’CfcJ  occurs.  Furthermore,  including  k  in  any  other  or  adding 

the  singleton  { k }  will  not  produce  a  partition  in  ik  ,  *1  k. 

respectively.  If  iQ  <  i^  then  {C^  .C^, .  • .  ,C^J  can  be  made  a  partition  of  r 

by  including  the  singleton  { k }  between  C.  and  C.  or  by  adding  k  to  one  of 

*0  1 1 


C.  ,C.  ...... C.  .  But,  we  have  seen  that  each  of  these  partitions  of  f 

1 0  *0  ll 

have  limit  0  as  w  -»  0,  except  for  the  two  that  include  k  either  with  C.  or 
k  *0 

C.  .  The  sum  of  the  probabilities  of  these  two  decompositions  approaches 
1  1 

P({C1«C2 . C? } )  as  wk  -*•  0.  Hence,  P(2.,k;w)  -*■  P( Jt,  ,k-1  ;  w  as  wk  -*■  0, 

and  the  proof  is  completed. 

•  We  now  return  to  the  consideration  of  the  sharpness  of  the  results  in 
Theorem  2.5.  More  will  be  established;  the  assumption  that  r  has  at  most 
one  minimal  or  at  most  one  maximal  element  will  be  relaxed.  However,  as  we 
have  seen,  the  bounds  given  there  are  not  sharp  if  r  is  decomposable.  We 
will  show  that  if  F  is  indecomposable,  then  the  bounds  riven  in  Theorem  2.5 
can  be  obtained  as  the  limit  of  P{SL,k;w^)  as  n  ■*  «°.  This  shows  that  the 

bounds  given  in  Theorem  2.5,  as  well  as  those  derived  from  them  by 
convolutions,  are  sharp.  Another  lemma  is  needed. 

Lemma  2.7.  If  (r,<<)  is  indecomposable  with  card.  (D  >  I,  then  there 
exists  a  j  e  f  with  f  -  ( j )  indecomposable  wrt  the  restriction  of  «  to  f  - 


element  in  B.,  then  <J>  y  (L-L.  )  y)  B  f  B.  implies  that  $  ^  (L-L  .JO  (B  - 

J  J  J  J  J  J 

{ k } )  ^  B^-{k)  and  hence, 

P(1,card.  (B.);w(B.))  ■*  P(1,card.  (B  {k} ) ;w(B fk}  ) )  as  w  ■»  0. 

J  J  J  J  K 

Thus,  P((BJ,BR,...,BjL})  -  P({B1  ,B2,...,Bj-{k} . B^} )  as  wk  -*•  0. 

Finally,  we  consider  decompositions  ino(^.  If  k  is  a  minimal,  but 
not  an  isolated  element  by  B ^ ,  then  <j>  j  (L-Lj_1  )  D  B^  implies  that  (L- 

Lj_1  )  0  ( Bj  -  { k } )  j  Bj-lk},  but  recall  L.  ^  U  {k }  £  £  .  For  ZR  <  0,  Zk//^  * 

-°°  and  so  B.  would  not  be  a  level  set  for  w  sufficiently  small.  Hence,  as 
J  K 

wk  -  0, 

P(1  .card.  (B.);w(B.))-P(Z  >0  )P  ( 1  .card  .  (B  {k}  )  ;w(B  {k} ) ) . 

J  J  *  J  J 

If  k  is  a  maximal,  but  not  an  isolated  element  of  B^ ,  then  a  similar 
argument  shows  that  as  wfc  *  0, 


P(1 .card.  (B  );w(B.)-P(Z  <0 )P( 1 .card . (B  {k| ) ;w(B . {k} )). 
J  J  K  J  J 


Because 


^1.2 . k-1 }  “  {L  ^  f1’2 . k-11:  L  ^},  every  partition  of 


T,  { B^  ,  , .  .  .  , B^ }  ,  can  be  made  a  partition  of  { 1  ,2 , . . , ,k-1 }  by  deleting  k 

from  the  B  containing  it.  Conversely,  let  (C  ,C0, . . . ,C„ )  y? .  ,  Set  i_ 

J  I  c-  £  X. ,  K- 1  U 

=  maxi  i :  there  is  an  a  e  C.  with  a  <<  k)  and  i,  =  minli:  there  exists  B  e  C  . 

]  1  l 


with  k  <<  B ) .  (Recall,  k  is  an  interior  element.)  Now,  i  <  i  ^ .  For,  if 


i,  <  iQ,  then  C,  0  ...  UC.^  . k.1).  whlG 


which  is  a  contradiction.  If 


iQ  =  i^,  then  [C^  fCR, . . . ,C^}  can  be  made  a  partition  of  f  by  including  k  in 
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decompositions  of  r  in  which  k  is  an  exterior  element  of  but  not  an 

(  \  ) 

isolated  element  in  B..  We  now  show  that  the  sum  in  (2.9)  over'**.. 

J  fck 

converges  to  zero  as  w,  -*  0.  There  are  two  cases  to  consider,  B.  =  { k j  and 

k  j 

card.(Bj)  Z  2.  If  «  { k } ,  then  because  k  is  an  interior  element  in  r,  1  < 
j  <  t.  As  w  -*•  0,  Var(Av(B  ))  and  Var(Av(B.))  are  fixed,  but  Var(Av(B.))  -» 

K  I  X,  J 

-.  Hence,  P(Av(Bt)  <  ...  <  Av ( B  ^ ) )  <  P{Av(B1>  <  Av(B^)  <  Av(B^))  -*•  0.  (See 

Robertson  and  Wright  (1982,  p.  305).)  This  also  shows  that  P(k,k;w)  +  0  as 

w,  •*  0.  We  next  consider  the  case  in  which  k  is  an  isolated  element  of  B. 
k  J 

and  card.  (B^)  -  2,  but  first  we  note  that  the  event  in  which  there  is  one 


level  set  in  the  projection  of  X(B.)  onto  the  cone  determined  by  <<_,  is 

J  bj 

determined  by 

Av(  (L-L.  .,  )  H  B.)  >  Av  ( B  . )  for  all  L  e&  with  <J>  f  (L-L  )  O  B.  J  B., 

J  J  0  J  J  J 

and  that  the  probability  of  equality  is  zero.  If  k  is  isolated  in  and 
card.  (B^)  >  2,  then  Lj_^  U  (k},  -  { k )  c aC  (if  a  <<  k,  then  a  e  L._^  or  a 

Q 

e  Bj ,  but  k  is  isolated  in  ;  and  if  k  <<  a,  then  a  e  L^^).  Thus,  P(1, 

card .  ( B .  ) ;  w( B  .  )  $ 

J  J 


P(X,  >  Av ( B .  )  and  Av(B.-{k})  >  Av(B.))  =  P(X,  >  Av(B.-fk))  >  X,  )  =  0. 
k  J  J  J  k  J  k 

Next,  we  consider  partitions  in3f^^.  Because,  B  -  { k }  j  <p ,  P(Av(B^) 

<  Av(B2)  <  ...  <  Av(B^) )  -  P(Av(B1)  <  ...  <  Av ( B j - { k } )  <  ...  <  Av(B^) )  as  wR 


*  0.  Replacing  Xk  by  Zk//w^,  it  is  easily  seen  that  if  k  is  an  interior 


19 


(k-1  ) 


with  P(l,k-1;w  )  the  probability,  under  «  ...  =  '  that 


— f  k-1 )  (k-1 ) 

u  has  exactly  8.  distinct  values.  (Of  course,  P(k,k-1;w  )  =  0.) 

Proof.  Assume  g  e  HQ.  As  we  have  seen,  the  MLSA  decomposes  T  into  disjoint 

sets  B.  ,B„ . B,_  ,  with  the  B.  successive  differences  of  a  nondecreasing 

1  2  h  l 

collection  of  lower  layers.  For  1  S  t  £  k,  let/ik  be  the  collection  of  all 
such  decompositions  with  h  =  £..  For  {B  ,B  ,...,B^}  c^>^k<  lnt  w(B  )  (X(B.)) 


be  the  vector  of  weights  w^  (means  X^.)  with  j  e  B^,  1  £  i  i  l.  Barlow  et 
al.  (1972,  eq  (3-23))  proved  that  P(SL,k;w)  = 

(2-9  . . V> 

with  P(  {  B1  .B^ . B^ } )  = 


P(  Av(B  )  <  Av  ( B  )  <  ...  <  Av  ( B. ) )  it  P(  1  ,  card .  ( B  .  )  ;w(  B .  ) ) . 

1  .  c  X.  11  • 

i=i 

Note  that  for  any  such  decomposition,  if  a  <<  8.  and  S  e  B.,  then  a  c  B^U  B?U 

...  vJ  B  ,  and  if  a  <<  B  and  a  e  B.  ,  then  6  e  B.  U  Further  note 

that  such  a  decomposition  does  not  occur  as  the  level  sets  with  positive 
l 

probability  if  it  P(  1  .card.  (B  )  ;v(B  ))  =  0. 

Let  i  denote  the  index  for  which  k  e  B..  Partition^.,  into  the 

J 

~J>  (  i  ) 

following  three  parts:  «•.  consists  of  those  decompositions  of  f  in  which 


-f  (?) 

k  is  an  isolated  element  of  consists  of  those  decompositions  of  T 

J  *k 


in  which  k  is  an  interior  clement  of  B.  and  £!,3)  consists  of  those 

J  ik 
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isotonic  wrt  «')CH|,  and  replacing  by  provides  a  lower  (upper) 
—2  —2  —2  —2 

bound  for  Xq1  and  and  E12^‘  As  we  have  Sfien  earlier  in  this 

section,  bounds  for  the  quasi  order,  <<*,  can  be  obtained  from  bounds  for 
its  reduction  <<**.  However,  <<’*  is  a  simple  tree  on  {1,2,...,e}.  By 
Theorem  1  of  Wright  and  Tran  (1985),  the  appropriate  bounds  for  this  simple 

0—1  —0+1 

tree  are  obtained  with  mixing  coefficients  1^)2  ,  which  is  our  A(it,e). 

The  proof  is  completed  by  applying  the  comment  given  after  the  proof  of  the 
Remark  in  this  section. 

If  T  =  T.  +  +  ...  +  r.  and  each  r  satisfies  the  hypotheses  of 

12  d  a 

Theorem  2.5  with  e  the  number  of  exterior  elements,  then  the  convolution  of 
a 

the  {A(t,e  )}  for  a  =  l,2,...,d  can  be  used  to  provide  bounds  analogous  to 

(2.8a,b,c,d) .  However,  this  convolution  is  Ad(£.,e)  with  e  *  e1  +  e„  +  ...  + 

e^ .  Now,  e  is  the  number  of  exterior  elements  in  r,  but  (A(St,e)  }  and 

{A^(t,e)J  are  not  the  same  unless  d  *  1.  In  fact,  if  d  >  1  and  c  >  0  the 

right  hand  side  of  (2.8  a,b)  ((2.8  c.d))  is  larger  (smaller)  if  A(S.,e)  is 
replaced  by  A^(l,e).  This  shows  that  if  one  first  decomposed  r,  then  better 

—2  —2  —2  —2 

lower  (upper)  bounds  for  xQ1  and  (x^  and  E^)  are  obtained. 

The-  proof  of  the  sharpness  of  these  bounds  rests  on  the  following. 

Lemma  2.6.  If  k  is  an  interior  element  in  (r,<<),  then 

lim  P ( l , k ; w )  =  P(l,k-1 ;w*K  ^ ) 

w,  -*0 
k 
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coefficients  which  provide  upper  (lower)  bounds  for  the  r^,  then  their 

convolutions  provide  upper  (lower)  bounds  for  r.  However,  if  the  breadth  of 

r  is  b  then  it  is  easy  to  show  that  b  =  b,  +  b„  +  ...  ♦  b.  is  the  breadth 
a  a  i  2  d 

of  T  and  the  convolution  of  the  (A.  (l,card(r  )-b  ))  is  { A  ( SL,k) )  .  Thus  it 

b  a  a  b 

a 

is  not  necessary  to  decompose  r  for  the  results  in  Theorem  2.2.  But,  as  we 
shall  see,  it  is  necessary  to  decompose  f  to  obtain  sharp  bounds  of  the  type 
provided  in  the  next  theorem.  It  should  be  noted  that  the  hypotheses  of 
Theorem  2.5  imply  that  T  is  indecomposable. 

Theorem  2.5.  Let  (f,<<)  be  a  partially  ordered  set  which  has  at  most  one 

maximal  element  or  at  most  one  minimal  element  and  let  e  denote  its  number 

of  exterior  elements.  For  w  e  HQ  and  any  w  with  w^  >  0  for  all  i  e  r, 

(2.8a)  PwCxo,  *  c]  2  AUie)P[Xf_1  *  c], 

(2.8b)  P  [E~.  >  c]  >  le.  AU,e)P[B  .  £  c], 

W  01  *=1  |(t-1),l(N-H) 

(2.8c)  2  c]  <  c)P[Xk_£  5  c],  and 

( 2 . 8d )  P  [E^n  >  c]  <  le0  A  ( l ,  e )  P  [  B  .  >  c]. 

W  i(k-l).l(N-k) 

Proof.  We  give  the  proof  for  the  case  in  which  T  has  at  most  one  minimal 
(.*  1  emont  .  The  proof  for  the  other  case  is  similar.  Because  T  is  finite,  it 
has  exactly  one  minimal  element.  We  suppose  that  the  minimal  element  is  1 
and  the  maximal  elements  are  2, 3.  ...,e.  Let  <<'  be  the  quasi  order  which 

requires  1  =  e*1  =  e+2  =  ...  =  k  <<*  j  for  j  =  2,3,-..,e.  If  x  e  R  is 

k 

isotonic  wrt.  <<’,  then  it  is  also  isotonic  wrt  <<•  Hence,  H]  *  |x  c  R  :  x 
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is  an  exterior  element  which  is  not  isolated.  Now  r  -  C  k }  with  the 
restriction  of  <<  has  breadth  b,  and  by  the  inductive  hypothesis  there  is  a 
weight  set  w  *  ( e ) )  for  which  |  P(  £  ,  k-1  ;  w '  (e  ) )  -  A^(£,k-1)J  <  e/2  for  2.  = 

1,2,...,k-1.  Set  w(e)  »  (w'fej.w^).  By  Lemma  2.3,  for  w^  small  enough, 

|P(t,k;w(e))  -  ^  P(K.,k-1  ;w*  (e) )  -  ~  P(£-1  ,k-1  ;w*  (e) )  |  <  e/2.  Using  the  fact 

that  Ab(2,k)  =  ^  A  (2-1,k-1)  +  ^  Ab(2,k-1),  we  see  that  | P ( 2 , k ; w ( e ) )  - 
Ab(2,k) | 

<  |  P(2,k;w(e) )~  PU.k-1  ;W(E))~  P(  2-1  ,k-1 ;  w*  ( e ) )  | 

♦  l|PU,k-1;w’(e))-A  U,k-1)|  ♦  i|P(2-1  ,k-1;w'(e))-A  (£-1  ,k-1  )  |  <  e. 

The  proof  is  completed. 

The  next  result  provides  lower  (upper)  stochastic  bounds  for  the  null 
—2  —2  —2  —2 

distribution  of  xQ1  and  E01  (x12  and  E12)  for  a  lar6e  class  of  partial 

orders  of  practical  importance.  This  class  includes  all  the  partial  orders 
discussed  in  the  Introduction. 

If  a  partially  ordered  set,  (r,<<),  can  be  partitioned  into  two 
nonempty,  disjoint  subsets  A  and  B  with  i  and  j  noncomparable  for  every  i  e 

A  and  every  j  e  B,  then  we  write  A  ♦  B  =  r  and  say  that  T  is  decomposable. 

Otherwise,  it  is  indecomposable .  If  r  can  be  written  as  the  union  of 

nonempty  disjoint,  subsets  r,  ,r. . with  the  elements  of  r  and 

1  2  d  a  B 

noncomparable  for  all  1  <  a  /  B  i  d  (ie.  r  +  r_  ♦  . . .  +  r  =  D,  then  the 

•  i  <1  a 

mixing  coefficients  for  r  are  the  convolution  of  those  for  r  and  r 

1  2  d 

(cf.  Barlow  et  al.  (1972,  p.  148)).  Hence,  if  one  obtains  the  mixing 


comparable  with  k.  (Such  exists  by  hypothesis.)  For  wk  sufficiently  small, 
Zk/«47^  >  max(X1 , . . .  ,Xk1  ) .  Using  the  MLSA  to  compute  the  projection  of 

(X  ,  Zk//wk),  we  see  that  the  first  i-l  level  sets  are  .B^, . . . ,B  ^ . 

Neglecting  the  set  with  probability  zero  on  which  Av(A)  *=  Av(B)  for  some  A 
and  B  distinct  subsets  of  r  and  noting  that  Av(A  (J  (kj)  -»  Av(A)  as  ♦  0 

for  ♦  f  A  CL  T  -  Ik),  we  see  that  for  sufficiently  small  wk>  L>  { k }  is  the 

next  level  set.  Hence,  for  Z,  >  0  and  such  w,  ,  M(X  ,  Z,  //w,  ;w)  = 

k  k  k  k 

-(k-1  )  (k-1  ) 

M(X  ;w  ).  The  proof  is  completed. 

Theorem  2 .  .  Let  (r,<<)  be  a  partially  ordered  set  with  breadth  b.  There 
exists  a  sequence  of  positive  weight  sets  w(n)  =  (w1 (n) , . . . »wk(n) )  for  which 

lim  P(g,,k;w(n))  =  A,  C 2. , k )  for  i.  =  1,2,...k. 
n-»"  b 

Proof.  The  proof  is  an  induction  on  k-b.  If  k-b  =  0,  then  the  partial 

» 

order  is  trivial,  there  are  no  comparable  elements.  Hence,  u  =  X. 

Suppose  that  k-b  >  0  and  by  relabelling  if  necessary,  let  (1,2,...,b) 
be  a  maximal  antichain.  We  show  that  for  an  arbitrary  c  >  0,  there  exists  a 
weight  set  w(e)  with  |P(C.,k;w(c) )  -  A^(i,k)|  <  e  for  t  »  1,2,.  ...k.  None  of 

the  elements  b+  1  , b*2, . . . , k  are  isolated  because  1 1 ,2, ....  b)  is  a  maximal 
antichain.  Furthermore,  we  will  show  that  at  least  one  of  b« 1 , h+ k  is 
an  exterior  element.  If  not,  then  they  are  all  interior  elements  and  one 
can  construct  a  chain  containing  b»l  with  both  its  smallest  ana  largest 
elements  in  {1,2,...,bJ.  But,  this  is  a  contradiction  since  no  two  elements 
in  {1,2,...,b)  are  comparable.  Without  loss  of  generality  w"  assume  that  k 


_  — ( k  —  1 ) 

projection  of  X  (X  )  onto  the  cone  determined  by  <<  (<<  restricted  to  F 


Lemma_2^_3.  Let  m1  -  v2  “  •••  “  !»•  If  k  is  an  exterior  element  but  not 


isolated,  then 


M(X;w)^  M(X(k_1  );w(k~1  *  I(  <b  p)(Xk) 


as  w,  ■»  0. 
k 


Proof.  Without  loss  of  generality,  assume  u  =  0.  We  suppose  k  is  a  minimal 
element.  The  proof  for  a  maximal  element  is  similar.  By  enlarging  the 

2 

probability  space  if  necessary,  we  can  obtain  Z  “  N(0,o_)  independent  of 

K  U 

_  _  _  £)  k_i  )  _ 

X1,...,Xk_1,  and  of  course,  M(X;w)  =  M(X  ,Zk//wk;w).  For  a  fixed 

element  in  the  underlying  probability  space,  we  consider  the  cases  Z.  <  0 

K 


and  Z  >  0.  If  Z  <  0,  then  for  sufficiently  small  w  ,  Z  //w~  < 

K  K  K  K  K 


min(Xj , . . . ,X  j  ).  Using  the  MLSA  and  the  fact  that  [kj  is  a  lower  layer,  we 

see  that  the  first  level  set  of  the  projection  of  (X  ,Zk//wk)  is  |k)  and 

— (k-1)  —( k -  1 ) 

the  others  are  the  level  sets  of  u  ,  the  projection  of  X  with 

(k-1  ) 

weights  w  .  Hence,  for  sufficiently  small  wk> 

M(X(k_1),  Zy/Suy-,v)  =  M(X(k_1  );w(k'1 })  ♦  1, 

Next,  we  consider  the  case  Zk  >  0.  Let  denote  tar-  level  sets 


-(k-1 ) 

for  u  and  let  B.  be  the  first  level  set.  containing  in  element 
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♦  b  is  a  stochastic  upper  bound  for  M  and  P[Y1  ♦  —  ♦  Y^  *  b  -  I]  - 

fk  ^  )2  kb  =  A. (£,k).  The  proof  is  completed, 
n-b  b 


We  next  show  that  the  bounds  in  Theorem  2.2  are  sharp.  Before  giving 
the  proof,  the  following  definitions  are  needed. 

Definitions.  An  element  i  e  T  is  said  to  be  maximal  (minimal)  if  there  does 
not  exist  j  e  r  with  i  <<  j  (j  <<  i).  Furthermore,  i  is  said  to  an  exterior 
element  if  it  is  a  maximal  or  minimal  element  and  otherwise,  it  is  an 
inter ior  element.  An  element  which  is  both  maximal  and  minimal  is  called 
isolated. 

The  minimum  lower  sets  algorithm  ( MLSA )  for  computing  projections  onto  a 
cone  of  isotonic  functions  is  used  in  the  proof  of  the  lemma  below.  A 
subset,  LC  T,  is  a  lower  layer  provided  i  e  T  implies  that  j  e  r  for  all  j 
<<  i.  Letc£  denote  the  collection  of  lower  layers  in  r.  For  ACT.  define 

Av  ( A )  =  L  .  w.X./L  ,  w.  and  set  L-.  =  0 .  Given  L„  CH  L.O  ...  C.  L  .  ,  the 
lc A  l  l  ic A  l  0  0  1  j 

MLSA  chooses  Lj  +  1  to  properly  contain  and  minimize  Av(L-L^)  over  all  such 


lower  layers  L,  and  in  fact. 


is  chosen  to  be  the  largest  such 


minimi zer . 


This  procedure  terminates  with  L^ 


r.  if  i  c  l.+1 


then  p.  =  Av(Bj).  The  are  the  level  sets  for  p  (cf.  Barlow  et  al.  (1972, 
p-  76)). 

—  (k-1)  —  —  —  (k-11 

Let  X  =  (  X  1  , X 2  ,  .  .  .  ,  X k _ 1  ) ,  w  =  (w1  ,w^ . wk-1^  and 


—  — ( k -  1  )  (k-1) 

M(X;w)  ( M  ( X  ;  w  )  )  be  the  number  of  distinct  values  in  the 
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H |  =  {x  c  R  :  x  is  isotonic  wrt  <<’}.  Examining  the  definitions  of 

2  2  2  2  o  — 

X01  *x\o>Eoi  and  Ei2  and  notin8  that  e12  is  an  incrcasinB  function  of  | { X  - 

—  2  —2  —2  —2  —2 
u||w/Q,  we  see  that  an  upper  (lower)  bound  for  Xq1  and  E^  (x12  and  E^)  is 

obtained  by  replacing  H1  by  Thus,  we  need  to  find  the  appropriate 

bounds  for  the  cone  Hj.  Robertson  and  Wright  (1982)  have  shown  that  because 

the  chi-squared  and  beta  variables  involved  in  (2.1)  and  (2.2)  ((2.3)  and 
(2.^))  are  stochastically  increasing  (decreasing)  with  ,  we  only  need  to 

find  DU.k)  with  £*  .  PU.k;w)  <  £*  DU.k)  for  j  -  1,2 . k.  If  M  is  the 

random  variable  which  denotes  the  number  of  distinct  elements  in  E  ( X | H *  )  , 

rk  —  2  rb  r 

we  seek  a  stochastic  upper  bound  for  M.  Since  i.  ^WjCX.  -  p^)  =  l  i c r 

a 

—  2  — 
w^  (X .  -  p^)  ,  we  may  obtain  the  projection  of  X  onto  H  ^  by  projecting  the 

subvector  of  X  with  coordinates  indexed  by  elements  of  r  onto  the  cone 

a 

which  requires  such  vectors  to  be  isotonic  wrt  the  restriction  of  <<’  to 

for  a  =  1,2,...,b.  Let  M  denote  the  number  of  distinct  elements  in  the  ath 

a 

such  projection  for  a  =  1,2,...,b.  Clearly,  M  are  independent 

and,  with  probability  one,  K  =  K  *  M  *  ...  4  M  .  Assure  p  t  H  .  Let  Y 

It  D  U  Ot 

B(card.(T  )  -  1,rr)  (B(0,i)  is  d«  ene  r  a  to  at  0)  for  1  ‘  u  £  b  with 

a  2  2 

Yj ,Y.}, . . . , independent.  Because  <<*  restricted  to  is  a  total  order,  Yq 

4  1  is  a  stochastic  upper  bound  for  y  for  a  =  1,2, ...,b.  Hence,  Y,  4  ...  4 

(1  I 


"small"  and  w 


•  •  •  •  f 


k-1 


are  not,  then  the  P(£,k;w)  will  be  close  to 


Ak_2(**k)  and  the  corresponding  bounds  would  provide  better  approximations 


than  the  equal-weights  ones-  On  the  other  hand,  if  w  ,  ....w^  are  small 
and  w^  and  wk  are  not,  then  P(£,k;w)  will  be  close  to  A(d,2).  (See  the 


proof  of  Theorem  2.8.) 

If  r  has  at  most  one  maximal  or  at  most  one  rinimal  element,  then, 
applying  Theorems  2.2  and  2.5,  we  see  that  for  c  >  0  and  u  e  HQ, 


Ye  (e  1  )2^  ep[v^  >  c]  <  P  [X2 

lx2.-i  _  J  ~  wLxoi 


rk  (k  b j  b  k  [  2 
H»b  i-b'  Lxi-1  J 


yk  k-b  b-k  (-  2  >  c]  <  p  >  c]  >  ye  (e-i)2i-ep[  *  >  c] 

Z-JL=bk JL-b ^  Lxk-£  J  wLX12  "  J  “  H=r 2.-1  ;  Lxk-H  “  J* 


e  ,e-1 .^1-e„r  2 


with  e  the  number  of  exterior  elements  in  f  and  b  the  breadth  of  r.  For  a 


simple  tree,  e  =  k  and  b  =  k-1;  for  a  simple  loop,  e  =  2  and  b  =  k-2;  for  a 

un'imodal  ordering  (ie.  1  <<  2  <<  ...  «  h  >>  h+1  >>  _ >>  k),  e  =  3  and  b  = 

2;  and  for  the  ordering  on  r  =  {(i,j):  1  <  i  <  R,  1  <  j  <  c)  given  by  (i,j) 

<<  (s,t)  provided  i  i  s  and  j  i  t ,  e  =  2  and  b  =  min(R.C). 

We  indicate  one  final  application.  Suppose  that  U,,U_,...,U  are 

1  2  m 


ordered  categorical  variables  of  interest  on  a  certain  population  and  that 
large  values  of  each  are  desirable.  Further,  suppose  that  one  wishes  to 
determine  whether  or  not  a  new  treatment  tends  to  produce  values  of 
which  are  at  least  as  large  as  under  a  control,  that  is  whether 


or  not  the  treatment  distribution  is  stochastically  larger  than  the  control. 

LRTs  for  the  case  m  -  1  are  discussed  in  Robertson  and  Wright  (1981).  Let 

p.  .  .  (q.  .  .  )  denote  the  proportion  of  the  population  with  U  at 

12  m  1  2  m 
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level  i^,  1  £  a  £  m ,  under  the  treatment  (control);  let  <<  be  the 

coord inatewise  ordering  on  r  ■  . i  ):  1  £  i  £  I  ,  1  <  a  <  m};  and 

i  2  m  a  a 

let  denote  the  lower  layers  wrt  <<.  Robertson  and  Wright  ( 1 9 7 ^  )  studied 
the  multivariate  concept  of  stochastic  ordering  which  requires 


H=  i,;  <  ,  ,  „  ,  P.-  5  *  I 


fo 


(i  ,i_,...,i  )  e  L  pi  i_...i  -  L{ i  , i», . . . , i  )  e  L  qi,i_...i 

12  m  12m  12  m  12m 

r  each  L  e  <£ .  H  is  equivalent  to  requiring  f(U,,U_ . U  )  to 

12  m 


be 


stochastically  larger  under  p  than  under  q  for  all  f:  R  -»  R  which  are 
nondecreasing  in  each  variable  with  the  others  fixed.  Let  A  denote  the  LR 
for  testing  H  versus  ~  H  based  on  independent  samples  and  T  =  -2  In  A.  It 
can  be  shown  that  the  significance  level  corresponding  to  t,  an  observed 
value  of  T,  ie.  sup  P  [T  2  t],  is  for  large  samples  approximately 

P»QGrl  p,C) 

(3.D  sup  I^Pd.kjpjPCx^  >  t]  with  k  =  1^...^. 

(The  details  will  be  given  elsewhere.)  The  breadth  of  this  r  is  b  = 
mindj  .I^..  •  •  »I  ),  and  applying  Theorem  2.2,  we  see  that  (3-1)  is 


yk  (k-bj  b  k  j-  2 

4e=b v  £.-b  ^  lx£-1 


tj. 
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